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We study the dynamical generation of masses for fundamental fermions in quenched quantum 
electrodynamics (qQED) at finite temperature in the bare vertex approximation, using Schwinger- 
Dyson equations (SDE). Motivated by perturbation theory, a further simplification is introduced by 
taking the wave function renormalization to be unity. In the zeroth mode approximation, the SDE 
for the fermion propagator resembles QED in 2+1 dimensions (QED3) at zero temperature with an 
effective dimensionful coupling a' = aT. For a fixed temperature, mass is dynamically generated 
above a certain critical value of this coupling. As expected, raising the temperature restores chiral 
symmetry and fermions become massless again. We also argue that by summing over the frequency 
modes and under suitable simplifications, qualitative aspects of the result do not undergo significant 
changes. 
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I. INTRODUCTION 

The perturbative treatment of gauge field theories at 
zero temperature is valid only for weak coupling, which 
fails to explain certain non-perturbative aspects of these 
theories such as dynamical mass generation for fermions. 
Schwinger-Dyson equations (SDEs) provide us with a 
natural platform to study this phenomenon because their 
derivation makes no assumption about the strength of the 
interaction. In quantum electrodynamics (QED) at zero 
temperature, it is well known that above a certain crit- 
ical value of the coupling, fermions can acquire masses 
through self interaction without the necessity of a non 
zero bare mass. This fact is in marked contrast to pertur- 
bation theory which prohibits the appearance of masses 
in this fashion as long as the bare mass is zero even if 
we did an all order resummation using renormalization 
group equations. In a thermal background, temperature 
is another parameter which comes into play and we ex- 
pect criticality in it too, triggering mass generation below 
a certain temperature and restoring chiral symmetry if 
we go above this value. The relevant SDE to investigate 
such a behavior is that of the fermion propagator. 

SDEs are an infinite tower of coupled equations relat- 
ing various Green's functions to each other. The two- 
point Green's function of the fermion propagator is cou- 
pled to all higher order Green's functions. Therefore, a 
non-perturbative truncation is essential to convert it into 
a solvable form. In quenched quantum electrodynamics 
(qQED), a favorite starting point is to make an ansatz 
for the fermion-photon vertex and then study the fermion 
propagator equation in its decoupled form. Inspired 
by some key features of gauge theories such as Ward- 
Green- Takasahi identities, Landau-Khalatnikov-Fradkin 
transformations, multiplicative renormalizability of the 
fermion propagator, etc., significant advances have been 
made in the hunt for a reliable vertex ansatz at zero tem- 



perature and the corresponding study of dynamical chi- 
ral symmetry breaking. However, for finite temperature, 
this venue has been much less explored, partially because 
of the richer tensor structure of the propagator and the 
vertex involved [1 H H E E The case of three di- 
mensional QED (QED3) at finite temperature has been 
studied in relatively more detail in the bare vertex ap- 
proximation [1 ©, El El El El El 

In this paper, we undertake the study of QED4 in 
the quenched approximation at finite temperature. We 
study the corresponding Schwinger-Dyson equation for 
the fermion propagator in the bare vertex approxima- 
tion. Under the additional assumptions of setting the 
wave function renormalization to unity and picking out 
only the most dominant frequency mode, the theory re- 
sembles QED3 at zero temperature (see also 15])- The 
effective coupling aT is dimensionful. In order to pursue 
the resemblance of the theory with QED3, we can take 
a' = aT along with T as independent parameters. The 
only notable difference is a temperature dependent mass 
term. The effect of this term is such that keeping the 
temperature constant, the increase of the coupling be- 
yond a certain critical value makes the fermion to acquire 
a mass. Conversely, increasing the temperature above a 
threshold value, we observe chiral symmetry restoration. 
We also confirm this behavior analytically. Moreover, in- 
creasing the value of the coupling, the critical value of the 
temperature also gets multiplied by the same factor. Wc 
must emphasize that the existence of criticality for the 
coupling and temperature is in contrast to the behavior 
of QED3 at zero temperature where no criticality exists: 
if masses are generated for one value of the coupling, 
they are generated for all other values of the coupling. 
Although not explicitly considered, the analysis can eas- 
ily be extended to the case in which a and T are taken 
as the independent parameters. 

The work is organized as follows: In Sec. [H] we set up 
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the SDE for the mass function in qQED with bare vertex 
at finite temperature using the imaginary-time formal- 
ism, starting from a simplified tensor structure for the 
fermion propagator. In Sec. IIIII we work under the as- 
sumption of considering only the contribution from the 
lowest Matsubara frequency finding a non-trivial numer- 
ical solution for the mass function. We study critical 
behavior of the mass function with the coupling strength 
as well as with temperature. With further simplifications 
in the integral equation, we also study this critical be- 
havior analytically, confirming the numerical results. In 
Sec. II VI we relax the assumption that considered the low- 
est Matsubara frequency as the single contributing mode 
and perform the sum over frequencies. We analytically 
study the resulting integral equation under the same sim- 
plifying assumptions as before and argue that, except for 
slightly different numerical coefficients, the behavior of 
the mass function appears to be qualitatively the same. 
We finally summarize and discuss our results in Sec. E] 



II. THE FERMION PROPAGATOR AND THE 
MASS FUNCTION 

The Minkowski space Feynman rules for the QED bare 
fermion propagator, photon propagator and fermion- 
photon vertex are, respectively, 



iS° F (P) = "4 



i 

J2 



Q 2 



(i) 



(hereafter we use capital letters to refer to four- vectors, 
whereas lower case letters are used for their components) 
where g is the electromagnetic coupling and £ is the usual 
covariant gauge parameter. Making use of the above 
rules, the SDE for the fermion propagator can be written 
as 



S F \P) = (S Q F )- 1 (P) 



«5 



d 4 K 



(27T) 



A r»(K,p)s F (K)>fn%,(Q),(2) 



where Q = K — P, T V {K,P) is the full fermion-photon 
vertex and Sf{P) is the full fermion propagator. 

At finite temperature, the most general tensor form for 
the fermion propagator can be obtained by noticing that, 
in addition to 1 and P , there are also two other Lorentz 
structures, namely p and P p. Therefore, we can write 



Sf(P) 



F 



P - bp - cPp - M 



(3) 



The Lorentz invariant functions F, b, c and M will in 
general depend on two Lorentz scalars P 2 and U ■ P, 



where is the four-velocity of the heat bath as seen 
from a general frame. We can choose these scalars to be 



Po = P^ 
p = [{P^U^f - P 2 ] 1 ' 2 . 



(4) 



The functions F and M are called the wave function 
renormalization and mass function, respectively. No- 
tice that for a parity conserving theory such as QED, 
c = O0. 

For our purposes, we will work in the approximation 
where 6 = and F = 1. The justification for this 
approximation stems from the fact that perturbatively 
F = 1 + 0(a) and b — 0(a). If a is small, we natu- 
rally expect F 1 and b w 0. Though it is well known 
that at zero temperature, inclusion of F is essential to 
recover the gauge invariance of the results such as the 
generated fermion mass [T^ |. it is also known that when 
ignoring F, the qualitative results, such as the shape of 
the fermion mass as a function of the coupling strength, 
are not changed. At finite temperature, we can expect 
that F and b play a similar role for the restoration of the 
gauge invariance 0. Nevertheless, to keep the discus- 
sion as simple as possible, here we will write the fermion 
propagator as 



S F (P) 



1 



P-M(P) ' 



(5) 



ignoring F and b in Eq. ©. 

Taking the trace of Eq. J5J, making use of Eqs. i ffl l 
and (J5J, and replacing the full vertex by its bare coun- 
terpart, the equation for the mass function can be written 
in Minkowski space as 



M{P) 



; (3 + 0c 

4-7T 3 



d A K 



M(K) 



I 



K 2 - M 2 {K) Q 



TT,- (6) 



We work explicitly in the imaginary-time formulation of 
thermal field theory, taking K — (iw„,k), P = (iu>,p) 
where u) n = (2n + 1)ttT and u> — (2m + 1)ttT are discrete 
fermion frequencies, with T being the temperature and 
n, m integers. In this formalism, Eq. © becomes 



M (P) 



(3 + fla , 

2tt 2 " 



M(K) 



[K-w) 2 + |k- 



k 2 + M 2 {K)} 



III. ZEROTH MODE APPROXIMATION 

A. The full equation 

A trivial solution to Eq. Q is M(P) = 0, which corre- 
sponds to the usual perturbative solution. However, we 
are interested in a non-trivial solution. To pursue it, as a 
first step we make the approximation where we take only 
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Mass Function in Landau Gauge 




T=0 - 
T=0.02 ■ 
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T=0.O35 
T=0.O37 - 
T=0.O375 - 
T=0.O379 



FIG. 1: The mass function M(p) for various values of tem- 
perature and a' = 1/47T in the Landau gauge. 



Generated Mass Versus Temperature in Landau Gauge 
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FIG. 2: The approximated generated mass M(0) for the 

fermion for various values of temperature and a' — 1/4-ir in 
the Landau gauge. 



the contribution of the zeroth mode, i.e., n = 0, in the 
sum over frequencies in Eq. Q and correspondingly take 
the same value for u>. We relax this approximation in the 
following section. It appears natural to take aT = a' 
and T as independent variables. After carrying out the 
angular integration, Eq. (0 becomes 



M(p) = 



(3+jK 
2ixp 



dk kM{k) 



TT 2 T 2 + k 2 +M 2 (k) 



In 



,(8) 



where p and k are the magnitudes of the vectors p and k, 
respectively. We solve Eq. (JHJ) by converting the integral 
equation into a set of simultaneous non- linear equations. 
The range of integration used is k = 1CP 3 to 10 3 . The 
number of points per decade to carry out the integra- 
tion was taken to be 31 which was found to be a suffi- 
ciently high number to obtain the accuracy with which 
the results have been quoted. In Fig. ^ we plot the mass 
function M for a' = I /Air and various values of T in the 
Landau gauge, i.e., £ = 0. With the increase of temper- 
ature, the generated mass decreases rapidly. The large 
p 2 behavior of the generated mass is depicted in Table I. 
We find that the mass function decreases as 1/p 2 in the 
asymptotic region of momentum. 

The generated mass can be approximated by M(0) 0, 
Il9| . since for small p, M(p) is rather constant. We draw 
it as a function of temperature in Fig. [21 in the Landau 
gauge. One can see a clear indication of the existence of 
a critical value T c of temperature above which there is no 
dynamical generation of mass. T c was found to be around 
0.037915. In order to see whether this critical value of 
the coupling can be obtained analytically, we analyze an 
approximated version of Eq. © in section IlTTCI 



B. Comparison with QED3 at zero temperature 

It is worth mentioning that at T = 0, Eq. © is the 
same as the one in QED3 at zero temperature in the 
approximations of using the bare counterpart of the full 
vertex and setting F = 1 (see Eq. (2.7) of Ref. 0), 
with the only difference that the factor 2 + £ has been 
replaced with 3 + £. The origin of this only difference 
is the the fact that the trace algebra is carried out in 
different dimensions in the two cases. We now mention 
an important contrast with QED3 at zero temperature 
where under identical assumptions, equation for the mass 
function reads as 



M(p) = 



(2 + 0" 



7Tp 



dk- 



kM(k) 
1 +M 2 (k) 



In 



k+p 



P 



(9) 



This equation is scale invariant. The solution for an ar- 
bitrary value of e 2 can be obtained from e 2 = 1 in the 
following fashion: 

M{p,e 2 ) = e 2 M(p/e 2 ,l) . 

Therefore, if chiral symmetry breaking solution exists for 
one value of coupling it exists for all values of coupling. 
There is no critical value of e 2 or a. A mathematically 
similar scaling law at finite temperature is 

M(p, e 2 , T) = e 2 M(p/e 2 , 1, T/e 2 ) . 

We have seen that for e 2 = 1, there is a critical temper- 
ature, T c = 0.037915, which separates the chirally asym- 
metric solution from the symmetric one, in the Landau 
gauge. Therefore, we have 



M(p, 1,21, 







V p 



for all temperatures T >c above its critical value. What 
happens if we double the value of coupling e 2 ? The scal- 
ing law implies that 



M(p,2,T) = 2M(p/2,l,T/2) 



4 



T p p 2 M(p 2 ) 


T P P 2 M{p 2 ) 


0.0 1000 3.89920 
640.4 3.89940 
304.699 3.89947 


0.037 1000 0.579892 
640.4 0.579921 
304.699 0.579931 


0.005 1000 3.84612 
640.4 3.84632 
304.699 3.84638 


0.0375 1000 0.386170 
640.4 0.386189 
304.699 0.386196 


0.01 1000 3.68662 
640.4 3.68681 
304.699 3.68687 


0.0376 1000 0.336004 
640.4 0.336021 
304.699 0.336027 


0.015 1000 3.41973 
640.4 3.41990 
304.699 3.41996 


0.0377 1000 0.277625 
640.4 0.277640 
304.699 0.277644 


0.02 1000 3.04294 
640.4 3.04310 
304.699 3.04315 


0.0378 1000 0.204152 
640.4 0.204163 
304.699 0.204166 


0.025 1000 2.54986 
640.4 2.54999 
304.699 2.55003 


0.0379 1000 0.0818793 
640.4 0.0818835 
304.699 0.0818849 


0.03 1000 1.92211 
640.4 1.92220 
304.699 1.92224 


0.03791 1000 0.0568280 
640.4 0.0568309 
304.699 0.0568319 


0.035 1000 1.08292 
640.4 1.08297 
304.699 1.08299 


0.037915 1000 0.0391043 
640.4 0.0391063 
304.699 0.0391070 



TABLE 1. The large p 2 behavior of M(p) for vari- 
ous values of temperature 



T p p 2 M(p 2 ) 


T p p 2 M{p 2 ) 


0.0 1000 2.41264 
640.4 2.41257 
304.699 2.41219 


0.028 1000 0.402201 
640.4 0.402190 
304.699 0.402126 


0.005 1000 2.35491 
640.4 2.35485 
304.699 2.35448 


0.0285 1000 0.268127 
640.4 0.268120 
304.699 0.268078 


0.01 1000 2.18138 
640.4 2.18132 
304.699 2.18098 


0.0286 1000 0.234034 
640.4 0.234028 
304.699 0.233991 


0.015 1000 1.89032 
640.4 1.89027 
304.699 1.88997 


0.0287 1000 0.194810 
640.4 0.194805 
304.699 0.194774 


0.02 1000 1.47509 
640.4 1.47505 
304.699 1.47482 


0.0288 1000 0.146417 
640.4 0.146413 
304.699 0.146390 


0.025 1000 0.905525 
640.4 0.905501 
304.699 0.905357 


0.0289 1000 0.0724110 
640.4 0.0724091 
304.699 0.0723976 


0.026 1000 0.762463 
640.4 0.762443 
304.699 0.762322 


0.02891 1000 0.0602952 
640.4 0.0602936 
304.699 0.0602841 


0.027 1000 0.600565 
640.4 0.600549 
304.699 0.600454 


0.02892 1000 0.0451881 
640.4 0.0451869 
304.699 0.0451798 



TABLE 2. The large p 2 behavior of M(p) for vari- 
ous values of temperature for the approximated log. 
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in Landau Gauge 
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FIG. 3: The approximated generated mass M(0) for the 
fermion for various values of temperature and two different 
values of coupling, a' = 1/47T and a' — 1/2tt in the Landau 
gauge. 



Let us substitute T — > 2T C in this identity which gives 

M(p,2,2T c ) = 2Af(p/2,l,T c ) . 

Thus M(p,2,2T >c ) = 0, i.e., the value of the critical 
coupling gets doubled by doubling the value of coupling 
e 2 . To confirm this scaling law we study Eq. (JSJ) for a' = 
1 /2ir and find that the new value of critical temperature 
is indeed roughly 0.075825. This is shown in Fig. 3. 

C. The equation with approximated log 

We can simplify Eq. fBl by employing the approxima- 
tion m, 



In 



y 



y 



= ^e{x-y) + —9{y-x), (10) 
x y 



which is a good approximation for x << y and x >> y. 
Moreover, for a small generated mass, we can also neglect 
the factor M 2 (k) in the denominator. Equation (JHJ) can 
then be written as 



M(p) = 



(3 + £K 



M(k) 



TT 2 T 2 + k 2 



9(k-p) 



+^e( P -k) . (ii) 



p 



If we repeat the numerical analysis of the previous sub- 
section, we find that the qualitative results do not un- 
dergo any drastic change. The mass function still lowers 
its height with temperature and above a critical value 
of the temperature (now 0.02893, see Fig. 5), only the 
trivial solution exists, i.e., M(p) = 0. This behavior is 
shown in Fig.0|in the Landau gauge. Moreover, the large 
momentum behavior of the mass function is again 1 /p 2 



as shown in Table II. This behavior for asymptotic values 
of momentum can be confirmed by converting Eq. i|ll|) 
into a second order ordinary differential equation 



d_ 

dp 



V 



;d M(p) 
dp 



4{3 + £)a'M(p) 







(12) 



with boundary conditions 
3 dM(p) 



P 



dp 







and 



M (p) = 



p=0 



p — >oo 



where we have dropped ttT in comparison with p. The 
solution of this equation is 



M(p) 



8(3 + 0"' 



Tip 



C1J2 4 



-c 2 Y 2 4 




(3+jV 

Tip 

(3 + £K 

Tip 



,(13) 



where J(x) and Y(x) are the Bessel functions of the first 
and second kind, respectively. The second boundary con- 
dition imposes c 2 = 0. As J(l/^/p) — ► 1/p when p — > oo, 
we conclude that M(p) — > 1/p 2 in this limit, as found nu- 
merically. However, such an analytic treatment washes 
out the information about the effect of temperature and 
we cannot find the critical value of coupling. To retain 
the effect of temperature, note that the factor irT in 
Eq. (|llfl serves as an infrared cutoff and therefore, we 
can rewrite this equation as 



M(p) = 



(3 + 



M(k) 



T ^ k 2 



8(k-p) 



k 



+ -6{p-k) .(14) 



P 



The corresponding differential equation is exactly as we 
had obtained before. Only the first boundary condition 
changes to 



;dM(p) 



dp 



= 



(15) 



This boundary condition leads us to the following equa- 
tion at the critical coupling 



where 



J 2 (4a) + a[Ji(4a)- J 3 (4a)] = 0, (16) 



= ^(3 + 0"'/^. (17) 



Note that a being a constant [determined by Eq. ((Ttjl) ]. 
we see that at the critical coupling, a' and T are pro- 
portional to each other, i.e., if we increase a' , the crit- 
ical value of coupling will increase by the same factor. 
The first non-trivial zero of Eq. (|l(jfl yields a = 0.958 
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FIG. 4: The mass function M(p) for various values of tem- 
perature and a = l/An for the approximated logarithm in the 
Landau gauge. 



Generated Mass Versus Temperature in Landau Gauge 
(With Approximated Log) 



-0.01 0.01 0.02 0.03 0.04 0.05 0.00 0.07 

Temperature 

FIG. 5: The approximated generated mass M(0) for the 
fermion for various values of temperature and a — 1/47T for 
the approximated logarithm in the Landau gauge. 



for a' = 1/4tt. This implies T = 0.02636 in the Landau 
gauge, which is within 10% of the actual numerical value. 

What happens when we are in a gauge other than 
the Landau gauge? In principle, the critical values a c 
and T c should be gauge independent because they trig- 
ger a phase transition. However, we do not expect it 
to be the case in the present work primarily because of 
the following reasons : (1) We assume the wave func- 
tion renormalization F(p) to be 1. The studies at zero 
temperature indicate that the inclusion of the equation 
for F(p) plays a crucial role in the partial restoration of 
gauge invariance, |18| . (2) The rainbow approximation 
neither respects the Ward-Green- Takahashi identity nor 
the Landau-Khalatnikov-Fradkin transformations, both 
being the consequences of gauge invariance. The gauge 
dependence of T c as a result of the numerical and ana- 
lytical analyses has been plotted in Fig. @. As another 
confirmation of the analytical result, Eq. I|17fl . we find 
that in our numerical calculation, the critical tempera- 
ture in an arbitrary covariant gauge, T c (£), is obtained 
from the one in the Landau gauge through the relation 
T c (0 = [(3 + 0/3]T c . 



IV. SUMMING OVER MATSUBARA 
FREQUENCIES 

Here we want to show that, under some simplifica- 
tions to analytically handle the expressions, the results 
obtained in the previous section appear to remain quali- 
tatively unchanged, when one considers the contribution 
from Matsubara frequencies other than the lowest one. 
Let us stress that this point is by no means evident since 
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Critical Temperature Versus Gauge Parameter Analytical 



Gauge Parameter 

FIG. 6: Critical temperature T c versus the gauge parameter 
£. The upper and lower curves correspond to numerical and 
analytical results, respectively. 



there are known examples in QED for which the contri- 
bution from the zeroth mode is exactly canceled by the 
contribution of the rest of the modes. One such example 
is provided by the one-loop fermion damping rate at finite 
temperature, which, as expected on physical grounds, 
when evaluated on the free mass shell vanishes due to the 
absence of phase space. This result can also be seen to 
arise as a consequence of the exact cancellation between 
the (divergent) contribution from the zeroth Matsubara 
mode and the rest of the modes [2l| . 

Since we will be primarily interested in order of mag- 
nitude estimates, we will carry out some rough approxi- 
mations aimed to preserve the possibility to analytically 
extract the behavior of the mass function in the critical 
region. The starting point is the sum over Matsubara 
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frequencies indicated in Eq. Q which we rewrite here as 
J(»w,p) = T^M(X) 

1 

X K + k2 + M2(A-)][K - u)2 + |k - p| 2 ] ■ 

(18) 



As a first approximation, we neglect a possible depen- 
dence of M(K) on the frequency, writing M{K) — M(k). 
In this approximation, the sum in Eq. I|18|l can be com- 
puted using standard methods [z| with the result 



I(iui, p) 



1 



4£ , i£ , 2 



1 + f(E 2 ) 



1 



i 



Ex + E 2 



f(E 2 )+f(E l ) 
1 

— i?i + i?2 «w + Ei — 



(19) 



where E x = y/k 2 + M 2 (k), E 2 = |k - p| and 

1 



/(*) = 
/» = 



1 

e*/ T + 1 



(20) 



are the Bose-Einstein and Fermi-Dirac distributions, re- 
spectively. Equation (|19fl contains the usual vacuum and 
matter contributions present at finite temperature, these 
last weighed with the statistical distributions correspond- 
ing to the boson and fermion internal lines in the self- 
energy diagram. 

Having carried out the sum over frequencies, it is 
now possible to perform the analytical continuation of 
the function I{iu>, p) to real energy values. We con- 
sider the analytical continuation to the retarded func- 
tion J r (po,p)- Although the function I r (po,p) contains 
in principle an imaginary as well as a real part, in order 
to make the connection with the analysis in Sec. 11111 here 
we concentrate on the real part. As stressed in Ref. Q, 
attention must be paid to the imaginary part when con- 
sidering the complete tensor structure of the fermion self- 
energy and/or a more complete form of the gauge boson 
propagator. The real part of the function I r is defined 
as 



Re/> ,p) = ~[/(^ 



Po + ie) + I (ho — > po - ie)] 



(21) 



and written explicitly as 



1 + f(E 2 ) - f(E 1 



x V 



1 



1 



p a - Ei - E 2 po + Ei + E 2 < 
f(E 2 ) + f(Ei) 



x V 



1 



Po - Ei + E 2 po + E 



?i-E 2 ) }' 



(22) 



where V denotes principal part. It is now easy to see that 
the integral in Eq. (JJJ can now be expressed in Minkowski 
space in terms of the functions 



h(Po,p) 



h(pa,p) 



d 3 k M(k) 
(2tt) 3 4EiE 2 
1 



pa + Ei + E 2 po - Ei - E 2 
f d*k M(k) r 
V J (2^4i^ /( M 



1 



1 



pa + Ei + E 2 po - Ei - E 2 
1 1 



h(po,P) 



pa - Ei + E 2 po + Ei - E 2 
f d 3 k M(k) - r 
V J j2^4Ei-E 2 fiEl) { 
1 



1 



pa - Ei - E 2 po + Ei + E 2 
1 1 



pa + Ei - E 2 po - Ei + E 2 



)} 



(23) 



which arc thus written to separate the different vacuum 
and matter contributions. The choice equivalent to work- 
ing with the lowest value of iu> is to take po = 0. Also, 
since we are interested in looking for critical behavior, we 
neglect M 2 for which E\ — ► k. 

As a second simplification, let us ignore the angular de- 
pendence of both the distribution functions or the mass 
function. This can be thought of as considering the sit- 
uation where p ~ 0. This is the region of interest for 
the present analysis since, as shown in the previous sec- 
tion, the small p region is the relevant domain for the 
development of a non-trivial solution for M. The above 
simplification allows to analytically perform the angular 
integrations. Therefore, keeping in mind that we are in- 
terested in the situation where p — > 0, the set of Eqs. ll'-il) 
becomes 



h(p) 
Hp) 
hip) 



1 



8tt 2 P Jo 
1 



p + 2k 



8tt 2 P jo 
1 



3tt 2 P Jo 



dkMik) In 
dkM(k)f(k) In 
dkM(k)f(k) In 



p- 2k 

p + 2k 



p-2k 
p + 2k 



(24) 



p- 2k 

Thus, in this limit, Eq. iQ can be written in Minkowski 
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space as 



CONCLUSIONS 



M(p) 



(3 + 0a 



2ttP jo 
1 + f(k) - />) 



dkM(k) 



In 



p + 2k 



p-2k 

We now use Eq. (jTHjl to write Eq. as 

(3 + Oar 



(25) 



M(p) 



2np 

p/2 



dfcM(fc) l + /(fe)-/(fe) 



4fc 
p 



+ / dfcM(fc) fl + /(*)-/(fc)l f). (26) 

./p/2 L J « J 

Deriving Eq. I|2tj|) and multiplying by p 3 , we obtain 
3 dM(p) 4(3 + 0" 



P 



dp 



p/2 



dfcM(fc) 1 + /(ft) - /(ft) ft. (27) 



For p — > 0, the integration region in Eq. (|27() is that 
of small k. Notice that, in spite of the presence of the 
Bose-Einstein distribution, the integrand is infrared safe 
since 



l + /(ft)-/(fe) 



k -> T 



(28) 



when fc — > 0. We can incorporate this behavior by 
considering the dominant contribution to the integrand, 
namely, the one coming from the low k expansion of the 
function /(ft) w T/k writing 



4 < 3+ ^t r /2 dk M(k). 



dp vi j T 

This requires the mass function M to satisfy 
,dM(p) 



P 



dp 



= 0. 



(29) 



(30) 



p=2T 



and, by further deriving Eq. J2HJ), to also satisfy the dif- 
ferential equation 



3 d 2 M(p) 2 dM(p) , 2(3 + £)a' 

- 6p 



dp 2 



dp 



M(p)=0. (31) 



Equation (|31|l becomes identical to Eq. I|12f) after the 
substitution T -> 2T. Thus, Eqs. and show 
that the scale for the onset of critical behavior can be 
obtained by considering only the lowest Matsubara mode 
in the analysis. 



We have studied the SDE for the fermion propagator 
to find non-trivial solution for the mass function in qQED 
under the assumption that the fermion-photon vertex is 
bare, the full fermion propagator can be approximated 
by Eq. (|SJ) and that only the zeroth frequency mode con- 
tributes dominantly. In order to stress the relation to the 
case of QED3, we have performed the analysis taking as 
independent parameters a' — aT and T. As expected, 
we observe that chiral symmetry is dynamically broken 
above a critical value of coupling a' for a fixed tempera- 
ture and is restored above a critical value of temperature 
while we hold the coupling constant. We also find that 
under simplifying assumptions, the result remains qual- 
itatively unchanged by summing over all the frequency 
modes. Although not explicitly worked out here, the 
analysis can be also made considering a and T as the 
independent parameters. 

In order to clarify the meaning of the mass function 
found in this work, recall that in vacuum, the ultravio- 
let behavior of the theory is connected to the dynamical 
symmetry breaking phenomenon. The connection, in the 
case of theories with anomalies (e.g. the chiral anomaly) 
comes, on general grounds, through the impossibility of 
simultaneously regularizing the theory while preserving 
the classical symmetries. In fact, it is well known that 
in qQED in 4 dimensions at T — 0, the generated mass 
function is proportional to the ultraviolet cut-off. How- 
ever, the important point to keep in mind is that this 
happens because this cut-off is the only mass scale avail- 
able in the theory. 

The former does not mean that the ultraviolet diver- 
gences of the theory are essential for dynamical sym- 
metry breaking. An example is provided by QED3 at 
T = 0, where chiral symmetry is dynamically broken 
despite the fact that the theory is ultraviolet safe and 
super-renormalizable. In this case, the mass scale avail- 
able is the dimensionful coupling. 

The above remarks apply to the case of the theory in 
vacuum and keep being true at finite T since the vacuum 
properties are left untouched by the temperature. 

We should stress however that the mass function found 
in this work is thermal, i.e. it vanishes for T = 0. This 
means that, unlike the vacuum case, the ultraviolet sector 
of the theory is not the one that provides the mass scale 
for dynamical symmetry breaking. This can be under- 
stood by noticing that any T = ultraviolet divergence 
is not present at finite T since the thermal distributions 
cut off the integrals at large momentum. Thus, the mass 
scale is provided by the temperature and in fact our re- 
sults show, for the case of qQED that we consider, that 
the generated mass function is proportional to T. 

Criticality arises due to the existence of a dimension- 
less parameter. In QED4 at T — 0, it is a. In contrast, 
in QED3 at T = 0, there is no such parameter and there- 
fore, chiral symmetry is dynamically broken for all values 
of the coupling. In QED4 at finite T, taking a' — aT to 
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be an independent parameter, the dimensionless parame- 
ter is provided by a' /T. Given the absence of ultraviolet 
divergences when considering only thermal effects, no ad- 
ditional scales come into play from the ultraviolet sector. 
The scaling law found is thus unaffected by the ultravio- 
let sector. 

It will be interesting to see the impact on results if 
we include the wave function renormalization F and the 
function b in which case the analysis should account for 
the imaginary part of the self-energy 0. It is also im- 
portant to repeat the calculation accounting for a more 
complete fermion-photon vertex, studying the impact of 
this treatment on the gauge dependence of the critical 
values for chiral symmetry restoration. All this is for 
future. 
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